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The standard derivations of electromagnetic energy and momentum in media take Maxwell's equations as the 
starting point. It is well known that for dispersive media this approach does not directly yield exact expressions 
for the energy and momentum densities. Although Maxwell's equations fully describe electromagnetic fields, 
the general approach to conserved quantities in field theory is not based on the field equations, but rather on the 
action. Here an action principle for macroscopic electromagnetism in dispersive, lossless media is used to derive 
the exact conserved energy-momentum tensor. The time-averaged energy density reduces to Brillouin's simple 
formula when the fields are monochromatic. The momentum density is not given by the familiar Minkowski 
expression D x B, even for time-averaged monochromatic fields. The results are unaffected by the debate over 
momentum balance in light-matter interactions. 
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I. INTRODUCTION 

Momentum balance in light-matter interactions is a subtle 
subject that has been debated for over a century. Discussions 
of the issues involved, together with references to much of the 
relevant literature, can be found in J2K6). This paper makes 
no attempt to contribute to the momentum-balance debate, in 
that there will be no investigation of the transfer of momen- 
tum (or energy) from light to matter or vice versa. Instead we 
isolate a related problem, but one that is mathematically well- 
posed and therefore capable of an exact and unambiguous so- 
lution. The problem solved here is: what is the conserved 
energy-momentum tensor of light propagating in a homoge- 
neous, dispersive, lossless magnetodielectric medium? This is 
a well-posed problem of macroscopic electromagnetism, but 
it appears not to have been fully addressed in the existing lit- 
erature on electromagnetic energy-momentum, where the pri- 
mary focus has been on momentum transfer between light and 
matter. Although the solution is an example of Noether's the- 
orem, the presence of dispersion means that the well-known 
expression for the energy-momentum tensor of a field theory 
in terms of its Lagrangian density is not valid in this case. 
In fact the result cannot be derived in any systematic way 
from Maxwell's equations or a Lagrangian density. Dispersive 
macroscopic electromagnetism must be formulated in terms 
of an action principle, and the Noether theorem for space- 
time translation symmetry must be derived from the action, 
leading to an exact conserved energy-momentum tensor. The 
problem in fact presents a significant calculational challenge, 
which may explain its absence from the momentum-balance 
literature. By isolating this problem from the conundrum of 
momentum transfer, we show that the issue of electromag- 
netic energy-momentum in dispersive media is amenable to 
some mathematically exact and unambiguous statements. 

The restriction to lossless media in this paper requires com- 
ment on the physical significance of the results. A complete 
treatment of macroscopic electromagnetism at all frequencies 



must include the absorption in the materials that is necessarily 
strong in some frequency ranges. Mathematically, this dissi- 
pation is a consequence of the restriction to retarded solutions 
of Maxwell's equations, which leads to the Kramers-Kronig 
relations HE). As is well known lf8l 4TOl . quantifying electro- 
magnetic energy-momentum in dissipative media is problem- 
atic, and in any case such energy-momentum will of course 
not be conserved. In optics, on the other hand, one often 
deals with limited frequency ranges where losses are negli- 
gible, but where material dispersion cannot be ignored. An 
example of the utmost experimental and practical importance 
is fiber optics ifTTl . where in many circumstances absorption 
can be ignored not only in the entire visible range, but also 
into the infrared and ultraviolet. In such a large frequency 
range dispersion plays a crucial role in light propagation and 
cannot be neglected IfTTl . The question of how much con- 
served energy-momentum is being transported by light has an 
unambiguous answer in these circumstances, and the result is 
experimentally significant. For example, a light beam encoun- 
tering an intense pulse can be frequency shifted through non- 
linearity of the medium, and the amount of electromagnetic 
energy that is frequency shifted depends, among other fac- 
tors, on the energy-momentum tensor of the light beam. One 
example of such an experiment is described in [12], where 
the light beam encountering the intense pulse was monochro- 
matic; this represents the simplest case, but for a light beam 
with a complicated spectrum the general expression for the 
energy-momentum tensor is required to fully describe the be- 
haviour of the beam. This exact energy-momentum tensor is 
derived here. 

The results obtained below generalize a well-known ap- 
proximation for electromagnetic energy in dispersive media. 
As is familiar from the textbook treatments HE), the Poynt- 
ing theorem that follows from Maxwell's equations does not 
directly lead to an expression for the electromagnetic en- 
ergy density in the case of lossless, dispersive media. The 
standard procedure (7[ [8] is to make a restriction to quasi- 
monochromatic fields, in which case one can extract an ap- 
proximate expression for the time-averaged energy density, 



due to Brillouin (Eq. (21 1 below). For general fields one is 
reduced to making a formal time integration of the Poynting 
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theorem [9|. But as noted above, in the dispersive case one 
cannot derive Noether's theorem from a systematic manipu- 
lation of the dynamical equations of the theory, as is done in 
the Poynting theorem; instead, one must start with the action 
that underlies the equations. Noether's theorem guarantees 
the existence of an exact energy density for arbitrary fields in 
lossless dispersive media and this energy density is derived 
here. 

In the case of electromagnetic momentum density, the dis- 
persionless result is the familiar Minkowski expression ||2]- 
6 1 D x B. We note again that this statement refers to 
the conserved electromagnetic energy-momentum tensor: the 
Minkowski momentum density D x B unambiguously gives 
the conserved electromagnetic momentum of light propagat- 
ing in a non-dispersive, lossless, homogeneous medium [2- 
6;j . It is shown below that in dispersive media the momentum 
density is more complicated, and is not given by Minkowski's 
formula. 



II. THE ACTION 

The results follow straightforwardly once an action princi- 
ple is written for electromagnetism in dispersive, lossless me- 
dia. For this purpose the dynamical variables must be taken to 
be the scalar potential cf) and vector potential A, defined by 
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-V(/>-d t A, B = VxA. 
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The relationship between fields in the time and frequency do- 
mains (the latter denoted by a tilde) is, for the example of the 
electric field, 
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where the reality of E(r, t) implies E(r, —ui) = E* (r, uj). In 
the frequency domain we have 



D = e £( r 7w)E, H = K K(r,o;)B, 



(3) 



where e(r, uj) is the relative permittivity of the (in general in- 
homogeneous) medium, kq = fj,^ , and the relative perme- 
ability is /u(r, uj) — k^.uj)^ 1 . As the medium can be as- 
sumed to be lossless in the frequency range of interest, e(r, ui) 
and k(t, uj) are real and even functions of uj O, and so they 
have the series expansions 

oo oo 

e(r, uj) = J2 £ 2„(r) uj 2n , k(v, u) = J2 «2n(r) " 3n , (4) 
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where these equations define the coefficients £2n(r) and 
K2n(r)- In practice, the series Q will represent a fit to the 
dispersion data of the material for the frequency range of in- 
terest; we take these series to be infinite, but the subsequent 
results also hold when they are finite series. Equations 



show that in the time domain the Maxwell equations (with no 
free charges or currents) are 

e V-[e(r s ift)E] = 0, (5) 
k V x [«(r,idt)B] = e e(r, id t )d t E, (6) 

while the other two Maxwell equations are automatically sat- 
isfied because of ([TJ. Note that series expansions of the form 
Q are standard in treating dispersion in frequency ranges 
where losses are negligible, for example in linear and non- 
linear fiber optics ifTTI . Moreover the transformation from the 
frequency to the time domain apparent in (|5]l-(j6]l, and the re- 
verse transformation, are a standard part of numerical solution 
procedures for the propagation of wave packets with spectra 
within the frequency range where expansions Q are valid, for 
example in the split-step method lfTTllT3i l. The action S[4>, A] 
for the potentials and A that gives the dynamical equations 
45)-© is 



S = 



r E-[ £ (r,£d t )E]-B-[/<r,id t )B] 



(7) 

Variation of cf> in |7]) gives the Maxwell equation (|5j, while 
variation of A gives ([6]). 

The fact that only even-order time derivatives occur in (|5]l- 
(|6]l is essential to being able to write an action principle for 
these equations. For dissipative media, terms with odd-order 
time derivatives occur in (|5]l-(j6]l, arising from the imaginary 
parts of e(r, uj) and n(r, ui) which are odd functions of ui 0; 
but the action |7]) would not generate these terms. In the ac- 
tion J7J, terms of the form E • <9 2n+1 E, for example, would 
not contribute to the dynamical equations because their varia- 
tion gives zero after integrations by parts: <5(E • (9 2n+1 E) = 
(SE) • <9 2n+1 E + E ■ <9 2n+1 <5E, and the second term becomes 
minus the first upon integrations by parts. It must of course 
be impossible to write an action in the dissipative case that is 
a functional only of <j> and A, since this would imply the ex- 
istence of a conserved electromagnetic energy in lossy media 
(see below). 



HI. ENERGY 

The action |7| is invariant under active time translations 
of the dynamical fields <f> and A, and this invariance im- 
plies, through Noether's theorem, the existence of a con- 
served quantity, the energy. The extraction of the conserva- 
tion law from the action is a standard technique of field the- 
ory (see lfl4l . for example). Even if one is unfamiliar with 
Noether's theorem, one can of course verify using the field 
equations |5]l-(|6]) that the resulting conservation law does in 
fact hold. The theorem [14| shows that if we make an ac- 
tive infinitesimal time translation <fi(r, t) — > cf)(r, t + £(r, t)), 
A(r, t) -> A(r, t + ((r, t)), but take the translation £(r, t) to 
vary in space and time, then the change in the action can be 
reduced to the form 



55= / d^xipdtC + S- VC) 



(8) 
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where p is the energy density and S is the energy flux, obeying 
the conservation law 



d t p + V • S = 0. 



(9) 



The variation 8S will clearly be linear in the infinitesimal 
translation £(r, t), but various numbers of integrations by 
parts are required to achieve the form ([8]l (surface terms pro- 
duced by the integrations are to be dropped). Use must be 
made of the following identities, which hold for arbitrary 
functions Y and Z, up to surface terms: 
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oo 2n 
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J2 E(-l) n+m £ 2 „(r)9r^ t 2 "- m Z5 t C 
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These identities were found by first working with small val- 
ues of n, where all the terms can be checked by hand, and then 



verifying the general expressions ( 10l-( 11 1 using Mathemat- 
ica. Identity (lOi with e(r,idt) replaced by n(r,idt) must 



also be used. In this way one attains the form <|8J> with 

Kq I 1 
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r(V^-fitA) • [e(r,id t )E] + B- [ K (r,id t )B] 
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S = - £ d t (j> e(r, £d t )E - K <9 t A x [/s(r, i<9 t )B] . 



(12) 
(13) 



It is straightforward to verify that ( |T2] > and ( [T3"] > obey the 
conservation law (|9|l when the fields satisfy Maxwell's equa- 
tions (|5]l-(|6]l. As in the case of dispersionless media, and in- 
deed vacuum, the energy density ( 12 1 and flux ( 13 i that di- 
rectly emerge from Noether's theorem are not gauge invari- 
ant ||7l fl4l . They are however equivalent to gauge-invariant 
quantities because they fail to be gauge invariant up to terms 
that identically satisfy the conservation law (|9j. Specifically, 
the quantities 



/V :=-s Q <t>e{r,id t )E i =:-f\\ 
/V := -« <Mr,^ t )(VM J ' - V* 'A*), 
identically satisfy 

9 t v l /V + v 4 (a f /V + v,/V) = o. 



(14) 
(15) 

(16) 



Comparing ( 16 1 with (pjj, we see that if V^/ 1 / is added to p, 
and dtPt + V 'j/V is added to S l , then the conservation law 



( 9 1 will still hold. Moreover, with use of Maxwell's equations 
( 5 i-Q, the energy density and flux that result from these ad- 
ditions are gauge invariant and are given by 



«0 p 

2 1 c 2 



p =— <; -E • [e(r, id t )E] + B • [«(r, t )B] 



1 oo zn 

- "a E E (-l)" +m £2n(r)9r _1 E • 9 t 2 "- m+1 E 

n— 1 m— 1 

oo 2n ^ 

+ ^ J] (-i) n+m ^n(r)ar _1 B • a 2 "- m+1 B }■ , 
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S =k E x [«(r,i9t)B]. 



(17) 
(18) 



The conservation law (|9|l for the final expressions ( fT7| and 
(TT~8b can also be verified using Maxwell's equations l5|l-(|6ll. 



can also be verified using Maxwell's equations ^5|-([6|. 
In (TT8 i we see that the energy flux is given by the Poynting 
vector E x H in the time domain, the obvious generalization of 



the non-dispersive result. The energy density ( 17 1, however, 
has no simple relation to the non-dispersive result, which only 
contains the first two terms in ( 17 1. It is a simple matter to 



express ( 17 1 and ( 18 i in the frequency domain. 



We can use the exact energy density (|T7j to derive the stan- 
dard textbook result [7 8] for the time-averaged energy den- 
sity of a monochromatic wave in a dispersive medium. Such 
a wave has an electric field 



E(r,<) = ^(E (r)e-^* 



c.c 



(19) 



and the B field is of the same form. When the monochromatic 



E and B fields are substituted into ( 17 1 and a time average is 
taken, all i-dependent terms vanish. The series in ( fT7| con- 
taining the E field is then 
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(20) 



4 ^ ^ --»\-v-u 4 

n— 1 m—1 n—1 



e de(r,uj ) 2 



The series in ( 17 1 containing the B field undergoes a similar 



simplification and the time-averaged monochromatic energy 
density p mono can be written 

e rf[wo£(r,w )] ,2 . p d[ui p(r,ujo)} . |2 

Anono - ~j H>0 + -j j -HO ) 

4 duio 4 dwo 

(21) 

which is Brillouin's formula Q O. For homogeneous 
media the monochromatic wave is a plane wave and it 
is easy to show, by dividing the time-averaged energy 
flow S by p, that the electromagnetic energy of the plane 
wave moves through the medium at the group velocity 
c/[d{u)Q -y 7 ' e{ujq)p{u)q)) I 'du>o\- It is of course expected that the 
energy of a monochromatic wave in a lossless, homogeneous 
medium should move at the group velocity, but here we have 
derived this fact for electromagnetism in lossless media with 
arbitrary dispersion. 
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permittivity is negligible; the real part is 






FIG. 1: The temporal profile of the electric field (top) and energy 
density (bottom) of a gaussian pulse in a medium with relative per- 
mittivity \22\ and fj, = 1. Note that the energy density is not zero 
at all the local minima; this feature is due solely to the complicated 
electric-field series in (17) . 



In the case of fields that are nearly monochromatic, Bril- 
louin's expression (21 1 can serve as an approximation of 
the time-averaged energy density, but only if e(r,ui) and 



/i(r, w) do not vary significantly over the range of frequency 
components in the fields |7|. As was recently pointed out 
in ifTUl . Brillouin's formula is also exact for time-averaged 
fields whose frequency components are uncorrelated, the im- 
portant example being thermal radiation. But the exact energy 
density for all fields in dispersive, lossless media is given by 

As mentioned in Section|II] the series representations (|4| on 
which our results are based will in practice be numerical fits to 
the measured dispersion data of the material in question, in the 
frequency range of interest. Clearly, a finite number of terms 
in the series will be sufficient for an accurate treatment of the 
dispersion, and the forgoing results hold in this case where all 
series terminate. As a theory exercise, however, we now con- 
sider an example of an infinite series expansion of the form (|4| 
for the permittivity, and verify that the expression ( 17 1 for the 
energy density is a well defined quantity. We take the stan- 
dard, damped harmonic oscillator model of a homogeneous 
permittivity Q. At frequencies well below the resonant fre- 
quency ujq of the oscillator, the imaginary part of the relative 



e(u) =1 



(uj 2 — ojq) 2 + 7 2 w 2 



n=0 



C-2 



n 



,2,,,2n 







m=0 



n + m 
2m 



2\ m 



(22) 
(23) 



The series expansion ( 22 1 converges for w < uiq and we con- 
sider a light beam in this material with a spectrum lying far 
enough below ujq for losses to be negligible. The evolution of 
the light beam in this dispersive medium can of course only be 
computed numerically, but we consider an initial, input gaus- 
sian pulse with central frequency ojq/2 and calculate its initial 
energy distribution. As is clear from ( 17 1, the effect of disper- 
sion on the initial energy density will be seen in the temporal 
profile of the pulse, rather than the spatial profile, so we con- 
sider a fixed point r = 0, which can be viewed as a boundary 
through which the pulse enters the medium. We take units 
with c = 1 and choose the values ojq = 10, SI = 30 and 7 = 1 



in the relative permittivity (22 1. The input electric field of the 
gaussian pulse centered on u = uiq/2 at the boundary r = 
is taken as 



E(*)|r=( 



cos(woi/2) 



(24) 



and is plotted in Fig. [T] We can numerically compute the in- 
put magnetic field B(<)| r= o by transforming to the frequency 
domain and use of the dispersion relation k = ^Jeiojc and 
the Maxwell equation iwB = V x E. To avoid inessential 
complications we ignore the transverse spatial profile of the 
pulse and consider a one-dimensional propagation (an actual 
one-dimensional propagation in a waveguide such as an op- 
tical fiber will involve an effective dispersion different from 
( 22 1, but our considerations here are purely for demonstration 
purposes). Using E(i)| r=u and B(£)| r=u we compute the in- 
put temporal profile p(t)\ T=0 of the energy density ( 17 i (with 
fj, = 1/k = 1) and the result is shown in Fig. [I] Note that no 
cycle-averaging has been performed; the exact energy density 
( 17 1 is not a cycle-averaged quantity. As well as the quanti- 
tative change in the energy density caused by the dispersion, 
there is a qualitative difference, visible in Fig.[T] compared to 
the same pulse in a non-dispersve medium. The local min- 
ima of the energy density in Fig. [T]occur at nodes of the elec- 
tric (and magnetic) field, but it is clearly seen that the energy 
density does not drop to zero at many of these local minima. 
This is purely an effect of the dispersion; in a non-dispersive 
medium these local minima are zeros of the energy density. In 
fact this effect of dispersion is due solely to the complicated 
electric -field series in (VT\ . The familiar E • D term in ( fTT] ) 
does not contribute to this feature, which can only be seen 



using the exact result ( 17 1. 
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IV. MOMENTUM 

Turning to the electromagnetic momentum, we note that 
momentum will be conserved only if the medium is homo- 
geneous, so that Noether's theorem applies to spatial trans- 
lations. It is instructive however to retain e and k as func- 
tions of position and thereby see how the momentum con- 
servation law fails to hold in the inhomogeneous case. To 
extract the conservation law associated with spatial transla- 
tion invariance we make the active infinitesimal translation 
0(r, t) -> cj)(r + T)(r, t),t), A(r, t) -> A(r + rj(r, t), t) in 
the action (|7]). Noether's theorem [14| shows that when the 
resulting change in the action is written in the form 



SS = - / d 4 x [Pidtrf +aPV j rf) , 



(25) 



then (in homogeneous media) Pi and a? obey the conserva- 
tion law 



(26) 



In this way we find the electromagnetic momentum density P 
and stress tensor a? . We will carry out this procedure with 
e and k varying in space; then ( |26| will fail to hold because 
extra terms will appear in this equation that contain spatial 
derivatives of e and k. Again, use must be made of the identi- 
ties (jTOf and ( fTTj ) to achieve the form ( |2"5j ) and the result is 



Pi = 



K 



r ViA j e(r,id t )Ej 



1 oo 2n 

-a E (-i)" +m £2n(r)ar- 1 £; i a t 2 ™- m v i ^ 



n—1 m— 1 

oo 2n 



^ (-l) n+ro «2n(r)flr~ 1 -B#-" , V < B J ' 



n—1 rn—1 



■X5/ +e Vi<p£(r,id t )Ei 
+ K V 1 A k K{v,id t )[V 1 A k - V k A j ], 



(27) 



(28) 



where C in (|2"8j ) is the Lagrangian density, i.e. the integrand in 
the action (|7j). To obtain gauge-invariant expressions for the 
momentum density and stress tensor, we note that the quanti- 
ties 



/V : ^A.dv.nhM-:' : / . . 



fV - 

/V :=-KoA-«(r,i9 t )(V^ fc 
identically satisfy 



* 3 
i 

V k A j ] 



(29) 
(30) 

(31) 



Thus, addition of Vjf* to P*, and of d t f\ j + V*/V' to f/. 
does not affect the momentum conservation law d26|. After 
these additions and use of Maxwell's equations (|5]l-(j6]), the 
momentum density and stress tensor are gauge invariant and 
are given by 

P = e e ijk [e(r,id t )&]B k 



oo 2n 

f EE 

n—1 rn—1 

oo 2n 

K 



{-l) n+m e 2n {v)dr- l E 3 dt 



2n—m 



e ^£(r, i$).EP - «o[«(r, id t )B,\B' 



1 



5/ [e 0J B fe£ (r,i9 t )E fc + K .B fcK (r, i9 t )B* 



(32) 



(33) 



where is the (completely antisymmetric) Levi-Civita ten- 
sor. It is straightforward to verify that, when the fields satisfy 
Maxwell's equations |5]l-([6]i, the momentum density (32 1 and 
stress tensor (33 i satisfy 



^•[Vi£(r,i5 t )]^' 



K 



B 3 [V^(r,za t )]B J , (34) 



so that the momentum conservation law ( p6*|) i ndeed holds 
for homogeneous media. The derivation of (|34|i in the non- 
dispersive case is familiar from the textbooks 1 15 1. 



Note that the electromagnetic stress tensor (33 1 is here de 



fined so that it is the spatial part of the electromagnetic energy- 
momentum tensor; a widespread convention for the stress ten- 
sor differs from this by a minus sign Q. The stress tensor 
( 33 1 is the obvious generalization to the dispersive case of the 
non-dispersive result; the fact that this expression describes 
the stress tensor in dispersive media is deduced in [8 | from 
completely different considerations. 

For the monochromatic wave ( fT9| ), the series in the mo- 
mentum density ( 32 1 containing the E field is, after a time 



averaging, 



oo 2n 

f EE 

n—1 rn—1 



ieinWa 1 1 ( 



EajViEg 3 



£ de(u} ) 
4 dojQ 



hn(E 0j ViE* j ). (35) 



There is a similar result for the series in (32) containing 
the magnetic field, and the complete result for the time- 
averaged monochromatic momentum density P in a homo- 
geneous medium is 



P =^eMRe(E x B*) - ^^Mlm^VES*) 



4 dujQ 



lm(B 0l VB* l ). 



(36) 



The Minkowski momentum density D x B would give only 
the first term in ( [36] ). With an additional restriction to a single 
plane wave of frequency ujq in the homogeneous medium, i.e. 



E (r) 



F c ik r 



Bo(r) 



1 

w 



kx E (r), 



w 



V / e(w )^(a;o), 



(37) 
(38) 
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FIG. 2: The temporal profile of the momentum density of a gaussian 
pulse in a dispersive medium. The pulse and dispersion are the same 
as those described at the end of Section [Til] and used for the energy 
density in Fig. [T] Note that the momentum density is not zero at 
local minima; this feature is due solely to the complicated electric- 
field series in i32\. 



the expression ( |36] l reduces to 



e(uj )n g (u} )k\£ \ , 



2c 

n g (u! ) = h p (ujo) +u> 



dn p (u! ) 

dojn 



(39) 
(40) 



where n 5 (wo) is the group index, n p (wo) = \/e(u}o)fi(uJo) 
is the phase index and k = k/|k|. Using the quantum pre- 
scription (for [a = 1) in lfl6l corresponds to taking |£o| 2 — > 
2Hlj / [£ n p ((jJo)n,g(uJo)] for a single photon; this gives from 
(39 1 a single photon momentum that is n p (wo) times the vac- 
uum value, in agreement with the canonical momentum in (6]. 
We note that the Minkowski expression (the first term in (f36*|) 



would give a correction factor of n^ujo) / n g {uJo) for the sin- 
gle photon [ 16 1, whereas the factor n p (uo) has been obtained 
here from the general momentum density (32). 



At the end of Section III we showed an example of the elec- 
tromagnetic energy density for a gaussian pulse in a disper- 



sive medium. For the same pulse and medium we use ( 32 1 to 
compute the input temporal profile of the momentum density, 



where the final series in the magnetic field in ( 32 1 is absent 
for this example. The spatial derivative of the electric field 



in ( 32 1 is computed by transforming to the frequency domain 
and using the dispersion relation k = y/eui/c (as discussed 



in Section III we assume one-dimensional propagation). The 
profile for the momentum density is shown in Fig. [2] and is 
seen to have the same shape as the energy density in Fig. [T] 
Similar to what was found for the energy density, the fact that 
the momentum density does not fall to zero at local minima is 
due solely to the complicated electric-field series i n p2] >; this 
feature can only be seen with the exact expression ( |32)> and is 
absent from the Minkowski contribution D x B in (|32|i. 

Electromagnetic angular momentum in dispersive media 
can be derived in a similar manner, though this will not be 
done here. The angular momentum density is not to be con- 
structed from the energy-momentum tensor but rather derived 
from the invariance of the action |7) under active rotations (for 
rotationally symmetric media). 



V. CONCLUSIONS 

We have derived the exact conserved energy-momentum 
tensor of light propagating in lossless, dispersive, homoge- 
neous media. Electromagnetic energy is also conserved when 
the medium is inhomogeneous and this case has been in- 
cluded. The energy flux and stress tensor in dispersive media 
have the same general form as in the non-dispersive case, with 
the permittivity and permeability becoming derivative opera- 
tors in the time domain. In contrast, the energy density and 
momentum density for arbitrary fields have no simple relation 
to the non-dispersive results. For time-averaged monochro- 
matic waves the Brillouin energy density is recovered, and the 



momentum density also takes a simple form (equation (36l). 
As stressed in the Introduction, the results do not address the 
question of momentum balance in light-matter interactions, 
which requires an analysis of momentum transfer between 
light and matter [2-6 1. The conserved energy-momentum ten- 
sor of light in lossless media is a well-defined quantity with its 
own experimental significance lfT2l : it has been derived here 
in complete generality within the framework of macroscopic 
electromagnetism. 
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